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ABSTRACT 

We present a detection of the Sunyaev-Zel'dovich (SZ) decrement associated with the Luminous 
Red Galaxy (LRG) sample of the Sloan Digital Sky Survey. The SZ data come from 148 GHz maps 
of the equatorial region made by the Atacama Cosmology Telescope (ACT). The LRG sample is 
divided by luminosity into four bins, and estimates for the central SZ temperature decrement are 
calculated through a stacking process. We detect and account for a bias of the SZ signal due to weak 
radio sources. We use numerical simulations to relate the observed decrement to I200 and clustering 
properties to relate the galaxy luminosity to halo mass. We also use a relation between brightest 
cluster galaxy luminosity and cluster mass based on stacked gravitational lensing measurements to 
estimate the characteristic halo masses. The masses are found to be around lO^^M©. 
Subject headings: cosmologyxosmic microwave background — cosmologyiobservations — galax- 
iesxlusters — Sunyaev-Zel'dovich Effect 
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1. INTRODUCTION 

New high resolution measurements of the sky at mil- 
hmeter wavelengths are enabling the use of the ther- 
mal Sunyaev-Zel'dovich (SZ) effect as a precise cosmo- 
logical tool. The effect arises from the inverse Comp- 
ton scattering of the cosmic microwave background 
(CMB) and high energy electrons, usually occurring 
in the hot intra-cluster medium (ICM) of gal axy clus- 
ters (ZeTdovich & Sunyaev 19691: ISunvaev fc Z ePdovichl 
ll97Qf ). For frequencies below 218 GHz, the distortion 
manifests itself as an arcminute-scale temperature decre- 
ment along the line of sight to a cluster. The amplitude 
of this decrement is nearly independent of redshift, mak- 
ing the SZ effect particularly useful for examining the 
high redshift universe. The effect can also be used to 
explore the role of baryonic physics in cluster evolution, 
as its amplitude is proportional to the thermal pressure 
of the ICM integrated along the line of sight. For a clus- 
ter in hydrostatic equilibrium, it is expected that the 
integrated SZ signal scales with total cluster mass. For a 
full descri ption of the SZ effe ct, see the r eview articles by 
iRephaelil (jl995f ) . iBirkinshawl (jl999f ) . and lCarlstrom et al.l 
(|2002f ). 

Analyses of large X-ray and optically selected clus- 
ter samples have recently illustrated the cosmological 
potential of cluster surveys (e.g., Vikhhnin et al. 2009; 
iMantz et all 120101 : IRozo et al.ll2010f ). The emergence of 
large- area SZ clu ster surveys, such as the South Pole 
Telescope (SPT; ICarlstrom et al.l I2Q09D and the Ata- 
cama Cosni ology Telescope (ACT; ' Fowler et all 120071 : 
iSwetz et al.l [2010). will supplement this previous work. 
Furthermore, upcoming results from the Planck satellite 
are expe cted to include a lar ge number of SZ-detected 
clusters (jBartlett et al.ll2008D . The SPT collaboration 
report ed its first SZ detections in iStaniszewski et aJI 
(|2009[ ) a nd has since identif i ed tw enty-two cluster candi- 
dates in 'Vanderlinde et al] (|2010l ) . Twenty-one of these 
clusters were optically confirmed in iHigh et all (^010). 
The X-ray properties of the SPT sample and the de- 
rived S Z /X-ray relations are presented in a follow-up 
study (jAndersson et al.l I2010D . The AC T collabora- 
tion p resented its first cluster detections in iHincks et al.l 
(120101) and has s ince reported on twenty-three clusters 
in Marriage et al. (2 011lV all of which have been opti- 
cally confirmed (Menanteau et al. 2010). The cosmo- 
logical implications of the ACT sample are discussed in 
ISehgal et al.l (|2QTTI) . 

In order to fully use cluster surveys for cosmological 
purposes, independent and robust estimates of the clus- 
ter masses are needed. Several studies have explored the 
relation between the integrated SZ signal Y and total 
cluster mass M for massive clu sters (M > 10^^M(7)) (i.e. 



Benson et aP 2004; Bonamen te et al.l [2QQ8I: iMelin et al l 



201 P; '"Huang et al.. i2010£ .Plagge et al 
grated SZ signal is defined as 



Wm . The inte- 



Y 



= / y{n)dn, 
Jn 



(1) 



where y is the usual Compton ^/-parameter and Q is the 
solid angle of the cluster. Assuming that thermal energy 
results solely from gravitational collapse, it is possible 
to derive self-similar scalings between observables (e.g., 
Y) and cluster mass. Furthermore, observations show 



that the SZ-mass relation is relatively insensitive to the 
specifics of cluster phys i cs (i.e., cooling, AG N feedback) 
(|Bonamente et al.ll2008l : 'Sehgal et al.l l2QllD . 

In this work, we report on the stacking of a subset of 
the ACT 2009 equatorial data at the positions of Lu- 
minous Red Galaxies (LRGs) measured by the Sloan 
Digitial Sky Survey (SDSS). Various studies have exam- 
ined the clustering properties of the S PSS LRGs (i.e., 
iZheng et al.]l2009l : iReid fc SDergelll2009D and have found 
that most of these galaxies reside in halos of typical mass 
M 10^^ - IO^^H-^Mq. This paper describes the stack- 
ing process and the detection of SZ signal after binning 
the sample by luminosity. Estimates for the character- 
istic halo mass of each luminosity bin are obtained us- 
ing two methods: gravitational lensing mass measure- 
ments and an em pirical model for halo bias as a func- 
tion of luminosity (Tegmark et al. 2004; Reid et al. 2010'; 
iReves e t al. 2008; Zehavi et al . 2005) that describe s the 
LRG sample well, as shown in fPercival et al.l (|2007[ ). Fi- 
nally, we present an initial analysis of the relation be- 
tween integrated SZ signal and halo mass. Unless oth- 
erwise stated, we assume a flat ACDM cosmology with 
Qrn 0.264, = 0.736, as = 0.80 and Ho = lOO/i km 
s-i Mpc-^ with h = 0.71 p^omatsu et al.l 12009). 

This paper is organized as follows. In Section [21 we 
describe the LRG sample and the ACT map. In Section 
[3l we describe the process of filtering, and in Section |4] we 
discuss our methods for binning and stacking the ACT 
data. We present the stacking results in Section [5] and 
examine the SZ-mass scaling relation in Section [6l We 
conclude in Section [71 

2. DATA 

In this section, we describe the subset of the ACT data 
used in this analysis (Sec. 12. ip and the LRG sample from 
the SDSS ([221). 

2.1. ACT Data 

ACT is a six-meter telescope operating in the Ata- 
cama Desert of Chile at an altitude of 5200 meters. The 
telescope site was chosen for its dry climate as well as 
the ability to observe in both the northern and south- 
ern skies. The telescope has three 1024-element ar- 
rays of transition edge sensors, one each operating at 
148 GHz, 218 GHz, and 277 GHz. For a more de- 
tailed introduction to the A CT instrument, observati ons, 
a nd d ata re duction, se e Fowler et all ( [200 7l ), Swetz et al I 
((2QT9 ) and iDas et all (|2QTT1 ). The science release from 
the ACT 2008 survey of a region south of the Galactic 
equator includes results on the CMB power spectruni 



and related parameter constraints (Fowler et al. I201C 
IDas et al. 2011; Dunkley et al. 2010; Hajian et al] l2010[ ). 
Results are also presented on compact millimeter sources 
(Marriage et al., 2010) and clusters (Marriage et alll2011l : 
[Menanteau et al.l [20101 : ISehgal et al.ll201lD . The present 
study uses data at 148 GHz from a subregion of the 
2009 equatorial survey. The subregion covers 268 square- 
degrees and lies between right ascensions 21^20"^ and 
03^40"" and declinations -01°30' and 01°30'. This region 
is coincident with the SDSS Stripe 82, which contains a 
rich multi- wavelength data set. Figure [His a map of the 
sensitivity across the subregion of the study, along with 
the locations of the LRGs. The map has been match- 
filtered for cluster detection as described in Section [3l 
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Fig. 1. — Sensitivity map with LRG locations. The map shows the sensitivity over the subset of the ACT 2009 148 GHz equatorial data 
considered for this study. The gray-scale encodes the noise rms in ^lK of a map match-filtered for cluster detection. The median noise in 
the map is 34 \iK in CMB temperature units. 

The median rms in the filtered map is 34 jllK in CMB 
temperature units. 



2.2. LRG Sample 

The LRG sample used in this study is from the SDSS, 
a survey that has p erformed five-band {ugriz) pho- 
tometry (jFukugita et a l. 1996) using a specialized wide- 
field camera ( Gu nn et al.l 11998 ) and multi-object spec- 
troscopy using a pair of spectrographs. The survey 
has imaged one quarter of the sky at high Galactic 
latitude and conducted spectroscopic follow up obser- 
vations of approximately one million of the detected 
objects ( Eisenste in etldl I2QQ1I : iRich ards et~aD "2002'; 
I Strauss e t al. 20021) • All data are processed automati- 
cally by pipelines that detect and nieasur e photometric 
properties of sources (|Lupton et a l.' '2001') . The SDSS 



has had seven ma jor data releases (A bazai i an et all 



2003, 



20081; 



[fooi, ^2005"; "Adelman- McCarthv et al.ii2006Tl2007r 
lAbazajian et al. 2009). 

The LRG catalog used in the curre nt work i s drawn 
from the sample prepared by Kazin et al.l (|2010f ) , which 
is taken from the seventh data release (D R7) of the SDSS. 
The d ata set was first developed by lEisenstein et al.l 
([2QQll ) and is publicly availableQ The DR7 LRG sam- 
ple contains ~1 10,000 galaxies and extends to redshift 
z 0.5. The subset used in this study contains 2681 
LRGs that fall within the 268 square-degree ACT map. 
The subset has a spectroscopic redshift range 0.16 < z < 
OAT and has rest frame ^-band ab solute magnitudes 
-23.2 < Mg < -21.2. K corrections f Blanton fc RoweisI 
^007) have been applied to the galaxies. The SDSS cat- 
alog of LRGs serves as a good tracer of large scale struc- 
ture, as LRGs are highly biased and trace more massive 
dark matter halos than other galaxy samples. They also 
exhibit a distinctive 4000 A break in their spectral en- 
ergy distributions, making redshift determination easier 
than for other galaxies. 

3. FILTERING THE ACT DATA 

To maximize the SZ signal of the clusters asso- 
ciated with the LRGs, we employ a matched fil- 

^ |http://cosmo.nyu.edu/^eak3Q6/SDSS-LRG.html| 



tei_jHaehnelt & Tegmari 119961 : iHerranz et all 120021 : 
iMelin et al. 2006). The filtering scheme used here closely 
follows Marriage et al. (2010), with exceptions noted be- 
low. At a position x in the map, we model the tempera- 
ture measured as the sum of the SZ signal plus the other 
components in the map (CMB fiuctuations, atmospheric 
fiuctuations, instrument noise and point sources): 



ST{x) = ^ 5To^ibi{x - Xi) + JTother(ic), 



(2) 



where STo^i is the peak amplitude of the i^^ LRG and 
bi is the unit-normalized ACT 148 GHz b eam function, 
taken to be isotropic with a FWHM of 1.4' (iHincks et al l 
1201 Of ). We first mask the brightest sources in the map, 
which amount to 146 objects. Only pixels associated 
with a source with SNR > 5 are masked. Thes e are pre- 
dominantly radio sources arriase eFaII[2QTQl ). None of 
these sources lies within 0.5' of an LRG position. Before 
applying the matched filter, we multiply the map, pixel- 
wise, by the number of observations per pixel normalized 
by the maximum number of observations in any single 
pixel, ^y A^obs {x) /A^obs,max • This Weighting accounts for 
local changes in the amplitude of the white noise and is 
equivalent to down-weighting pixels with a high white 
noise rms. This weight scheme is also used during the 
stacking process when determining the matched filter 
decrement STo (see Section |4?2]) . Next, the map is fil- 
tered in Fourier space using a matched filter: 



STmtik) 



b*{k) |(5rother(fc)|-^ 5T{k) 

Jb*{k') |5Tother(fcO|-' kk')dk'' 



(3) 



Since the power spectra of the atmosphere and the CMB 
scale as for / > 1000, the filter is approximately a 
fourth derivative operator convolved with the beam. The 
choice of the beam function b for the spat ial profile of the 
cluste r differs from the method used in iMarriage et al.l 
(|2011[ ), which employs a /3- model for the cluster profile. 
From comparisons with simulations, we expect the ACT 
148 GHz beam function to be a good match for the spa- 
tial profile of clusters and hence adopt it in the matched 
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TABLE 1 

Luminosity Bins and Halo Masses for Stacked LRGs 



Bin A^bin {Lo.ir) Lo.ir Range {z) M200/i 

1 26 30.1 25.9 - 44.9 0.39 4.26 ± 1.06 2.59 ± 0.87 

2 60 23.6 21.8 - 25.3 0.40 2.89 ± 0.20 2.05 ± 0.43 

3 674 16.7 13.8 - 21.8 0.37 1.71 ± 0.34 1.43 ± 0.31 

4 1921 10.0 1.1 - 13.8 0.30 0.80 ± 0.24 0.82 ± 0.22 



^Determined from analysis of halo bias. 

^Determined from lensing measurements in lReves et al.l (|2008l ) . 

filter. It is particularly effective at reducing scatter from 
CMB noise, to which broader filters are more sensitive. 

Rather than explicitly modeling the noise contribution 
^^^other(3^), we usc the original input map with sources 
masked for the power spectrum of ^Tother(^) in equation 
[3l As the model for STother{x) is not well- conditioned, 

we smooth |^Tother(^)|~^ with a Gaussian filter {a = 
5 Fourier space pixels). We set outlying Fourier space 
pixels with amplitude greater than 10 times the median 
pixel value to the median value. We also apply an addi- 
tional low-^ filter in Fourier space that grows from zero 
at ^ = to unity at ^ = 1200 as sin^(7r^/2400). This filter 
accounts for any modification of CMB statistics due to 
multiplying by the non-uniform weight map and makes 
our result more robust to any potential non-Gaussianity 
in the large-scale noise that might arise from residual at- 
mospheric noise in the maps. After finishing the filtering 
process, we do not include any sources that fall within 
10' of the map edge in order to mitigate noise from edge 
effects. This reduces the map size from 284 square de- 
grees to 268 square-degrees. The entirety of the LRG 
sample discussed previously (2681 objects) falls within 
the bounds of the reduced map, leading to roughly 10 
LRGs per square-degree. 

4. STACKING METHODS 

While ACT can easily detect massive clusters (M > 
lO^^Mo) (jMarriage et al.l 120 111 ). CMB fluctuations and 
detector noise limit the direct detection of less massive 
clusters. Since the mean of both the CMB fluctuations 
and atmospheric noise should be zero, we stack submaps 
around LRGs to enable the detection of weaker SZ sig- 
nals. 

4.1. Luminosity Bins 

Before stacking the LRGs, we bin the sample by r- 
band luminosity K corrected to 2; = 0.10, i^o.ir- Table [1] 
shows the four luminosity bins used in this analysis. This 
table also lists mass estimates for each bin, as described 
in Section [5^31 The luminosities have been computed us- 
ing: the Galactic e xtinction-corrected (using maps from 
iSchlegel et al.l ll998l) Petrosian magnitudes from the DR7 
CAS and their measured errors. They have been K cor- 
rected to z = 0.10 using; the kc orrect-v4_2 software 
package (jBlanton fc RoweisI 120071 ). There have been no 
added evolutionary corrections. We also refer to the K 
corrected ^-^r-band a bsolute m agnitude of an LRG as 
Mo.ir. Following Blan ton et a l. (2003), 

Mo.lr = mr - DM{Z, Qrn, ^A, h) - i^o.lrr(^), (4) 



where DM(z, ^7^, (^a, /i) is the distance modulus in a 
{Qrn = 0.3, Oa = 0.7, h = 0.71) cosmology and i^o.irr(^) 
is the K correction from the r-band of a galaxy at 
redshift z to the ^- ^r-band. Luminosities are reported 
in sola r luminosities, w here we have used M0^o.ir = 
4.76 (Blanton et al.| [2QQ3l) . These definitions of luminos- 
ity and magnitude will facilitate the estimation of halo 
masses for our binned data in Section 15.31 

4.2. Algorithm 

We stack filtered submaps with area 10' x 10' centered 
on each LRG in a given luminosity bin, and then esti- 
mate the matched filter decrement STq from the central 
temperature value in the stacked map. Recall that the fil- 
tered map is weighted by the number of observations per 
pixel normalized by the maximum number of observa- 
tions in any single pixel. To minimize th e errors induced 
by th e 0.5' pixel size of the ACT maps (Marri age et al.l 
I2010f ). the pixel size of the submap centered on each LRG 
is decreased using Fourier interpolation. For the present 
analysis, the pixel size is decreased to 0.03125' and the 
submaps are recentered. Before determining JTo, the 
stacked map for a given bin is convolved with the ACT 
148 GHz beam function b to recover decrement informa- 
tion from the surrounding pixels. Lastly, we sum over 
the five pixels within 4" of the centroid of the stacked 
map and divide by the sum of the weights to compute 
STq. We compute the weight for each submap in the 
stack before decreasing the pixel size and use the value 
of ^y A/obs/^obs,max at the Central 0.5' pixel. 

5. STACKING RESULTS 
5.1. Matched Filter Decrement STq 

Tables [2] and [3] contain estimates for the matched filter 
decrement 6T0 obtained for each of the four luminosity 
bins. Due to radio contamination (see discussion below), 
we report results for the full LRG sample as well as a 
radio-quiet subset, which excludes approximately 10% of 
the LRG sample. The value for dTo reported in Tables 
[2] and [3] is the temperature decrement in the ACT map, 
recovered from the matched filter. Figure [2] shows 7' x 7' 
stacked submaps for each luminosity bin for both LRG 
samples. Stacked maps obtained when using random 
LRG positions are also shown, and it is apparent that 
there is no clear detection in these maps, as expected. 

In addition to the central peak in the stacked images, 
several other significant peaks are evident, particularly 
the deep cold spot at the upper edge of the stacked map 
for Bin 2, which has a depth of around 24 jllK and is 
significantly deeper than the stacked cluster at the cen- 
ter. For this stacked map, the pixel noise is around 5 
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TABLE 2 

Central SZ Temperature Decrement and Integrated SZ Signal for All Stacked LRGs, Binned by Lo.ir 



Bin 


A^bin 


5To 


SNR 


^200p 




Y2m-p d\ E{z)-y^ 






laK 




arcmin^ 




Mpc^ 


1 


26 


-27.5 ± 11.3 


2.4 


(2.0 ± 0.7) xlO" 


-4 


(1.4 ± 0.6) xlO-^ 


2 


60 


-9.9 ± 5.3 


1.9 


(6.2 ± 2.7) xlO" 


-5 


(4.4 ± 2.4) xlO-^ 


3 


674 


-4.7 ± 1.5 


3.1 


(2.4 ± 0.7) xlO" 


-5 


(1.7 ± 0.6) xlO-6 


4 


1921 


-1.6 ± 1.0 


1.6 


(4.5 ± 4.3) xlO" 


-6 


(4.6 ± 2.8) xlO-'^ 








TABLE 3 






lATURE Decrement and Integrated SZ Signal for Stacked Radio-quiet LI 


Bin 


A^bin 


5To 


SNR 


^200p 




Y200P d\ E{z)-''r^ 






laK 




arcmin^ 




Mpc^ 


1 


21 


-28.2 ± 12.9 


2.2 


(2.2 ± 0.8) xlO- 


-4 


(1.5 ± 0.7) xlO-^ 


2 


51 


-10.6 ± 4.9 


2.2 


(5.4 ± 2.4) xlO- 


-5 


(4.7 ± 2.2) xlO-^ 


3 


587 


-6.1 ± 1.6 


3.8 


(3.3 ± 0.7) xlO- 


-5 


(2.3 ± 0.6) xlO-6 


4 


1732 


-2.1 ± 1.0 


2.1 


(1.1 ± 0.9) xlO- 


-5 


(7.9 ± 6.2) xlO-'^ 



Note. — The estimates for are determined from 5To using microwave sky simulations (Sehgal et aLII2010r ). They correspond to 

the halo masses M200P derived from considering luminosity-dependent halo bias (Table The uncertainty in 5To is determined from 
bootstrap resampling. The central decrements and associated errors are measured in CMB temperature units. 



|LtK, making the upper-edge cold spot a 5cr discrepancy. 
Examination of the 60 individual maps in the stack re- 
veals that two maps contain decrements of < 150 |llK at 
the location of the cold spot. In addition, 11 additional 
map sections show decrements larger than 50 |llK at the 
approximate position of the upper-edge cold spot. This 
large number of decrements in similar positions relative 
to the LRGs in our catalog clearly is not a random phe- 
nomenon, but rather reflects the fact that galaxy clusters 
and groups are strongly clustered on the sky. 

The errors in 5To presented in Tables [2] and [3] are com- 
puted for each bin using bootstrap resampling^ These 
errors measure the scatter in the properties of individual 
LRGs. Figure [3] shows the resampled distribution of the 
mean of 5To for each of the four luminosity bins for the 
full and radio-quiet LRG samples. The distribution of 
each bin is well- approximated by a Gaussian function. 

Table m contains the matched filter decrements 5To ob- 
tained when stacking various subsets of the LRG sample. 
These subsets are designed to test various systematic as- 
pects of the ACT data. The first four tests in Table H 

TABLE 4 

Matched Filter Decrement 5To for Various Systematic 
Tests of the ACT Data 



Test 


A^bin 


5To (laK) 


Top 


1295 


-2.51 ± 1.16 


Bottom 


1386 


-3.13 ± 1.10 


Left 


1406 


-2.94 ± 0.99 


Right 


1275 


-2.67 ± 1.35 


High A^obs 


1455 


-2.15 ± 0.98 


Low A^obs 


1226 


-3.85 ± 1.36 


Radio-loud 


290 


3.27 ± 2.65 


Radio-quiet 


2391 


-3.56 ± 0.84 



Note. — The significant negative signal shown is expected. The 
central decrements and associated errors are measured in CMB 
temperature units. 



^ We create 10,000 bootstrap realizations from the values of 5To 
in each bin and use the distribution of the resulting bin averages 
to estimate the uncertainty 



probe spatial discrepancies in the data by splitting the 
ACT map by right ascension and declination. We con- 
clude that the signal is consistent throughout the map. 
We also examine variations in white noise amplitude by 
splitting the sample spatially by the number of observa- 
tions A/obs per pixel. We determine the median value of 
Aobs per pixel and then stack subsets of LRGs that lie 
above and below this median value. We conclude that 
the signal is present for both high and low noise bins in 
the map. 

Radio source contamination at 148 GHz appears to 
have a significant effect. Using the SDSS DR7 CAS, we 
identified LRGs with radio co unterparts within 2 " in the 
I A GHz VLA FIRST Survey (jBecker et al.l[T995h and de- 
fine these LRGs to be radio-loud. There are 290 radio- 
loud LRGs in the sample: 5/26 in Bin 1, 9/60 in Bin 
2, 87/674 in Bin 3 and 189/1921 in Bin 4. Contamina- 
tion by FIRST radio sources is potentially a concern for 
~10% of the LRG sample. We also compute the matched 
filter decrement for each bin excluding the 290 radio-loud 
LRGs. While the stacked radio-quiet LRGs show a sig- 
nificant (> 4cr) decrement, the radio- loud LRGs show 
clear evidence of radio source contamination. For the 
rest of this analysis, we report results for both the full 
and radio-quiet samples (see Tables [2] and [3]) . 

5.2. Integrated SZ Signal Y 

We use numerical simulations ('Sehgal et al.l [2QTqI ) to 

relate the matched filter decrement 5To to the integrated 
SZ signal Y (equation [1]) . In an unfiltered map, where 
the cluster is larger than the beam, Y STqO'^^ where Oc 
is the angular size of the cluster. For a cluster of fixed 
physical size, Y ex STod'^^ where dA is the angular diam- 
eter distance to the cluster. Because the filter acts as a 
operator and removes large-scale signals in the map, 
the expected relationship for nearby clusters of fixed size 
is likely even steeper. However, rather than rely on this 
approximate scaling to weight the observations, we apply 
the same filter used on the ACT data to the simulated 
maps and compute 6T0 for each of the clusters. The sim- 
ulated maps are first convolved with the ACT 148 GHz 
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beam function. Then, motivated by the angular diameter 
distance scahng, we fit a model for the Y — STq relation 
which allows for a power-law dependence on redshift, 



TABLE 5 

Best-Fit Parameters for the Y2oop - STo relation, obtained 

FROM SIMULATION AND FITTED USING EQUATION [S] 



C ST,, 



(5) 



where Zi is the cluster redshift and Yest,i serves as an es- 
timate for the SZ signal within a disk of radius R200 for 
the i^^ cluster. We define R200 as the radius within which 
the average density is equal to 200 times the mean matter 
density of the universe, p{z) = pcrit{z = 0)flrn{^ + z)^- 
The redshift dependence in equation [5] for the i^^ cluster 
is determined mainly from the angular size of the cluster 
profile, after convolution with the matched filter. For a 
fixed value of JTo, lower redshift clusters appear larger 
on the sky and have larger values of Y than clusters at 
higher redshift, as shown in Figure |H As we are con- 
sidering clusters of comparable size to the ACT beam 
function (see Figure [2]), we expect the integrated SZ sig- 
nal Y to depend linearly on the matched filter decrement 

The relationship between Y and 6T0 is obtained for the 
mass range of each luminosity bin. Here the mass range 
is represented by the la uncertainty in the central M200P 
value, listed in Table [TJ Due to systematic uncertainties 
in determining halo mass from LRG luminosity, we per- 
form this analysis for two separate groups of mass bins 
(see Section The best fit values and corresponding 
errors for C and /3 are summarized in Table O The un- 
certainties reported are the formal la errors computed 
using Levenberg-Marquardt least squares minimization. 
The correlation between Yest and I200P is shown in Figure 
[5l The scatter in this correlation results from a combi- 
nation of noise in the simulated maps and the fact that a 
cluster's profile and the center of SZ signal is in general 
not that of the matched filter profile. 

The values for I200P are found for the j*^ luminosity 
bin by summing over each LRG in the bin as follows. 



Y 



200p,j 



(6) 



where 6To^i is the matched filter decrement recovered 
from the filtered ACT map and Nj is the number of 
LRGs in the j*^ luminosity bin. The weight Wi is related 
to the number of observations per pixel (see Section 
Values for I200P are reported in Tables [2] and [3l These 
values correspond to the halo masses M200P derived from 
considering luminosity-dependent halo bias. The error in 
^200p represents the uncertainty after propagating errors 
in ^To, /3, and C. These formal errors do not include 
the systematic uncertainties in our model w hich assumes 
that the gas profiles in Sehgal et al." ("2010") are accurate 
representations of the profiles for typical clusters. 

5.3. Estimates for Halo Mass M200P 

We take two different approaches to estimating the 
characteristic mass associated with halos of a given lu- 
minosity. The first approach is to use measurements of 
the halo bias, defined as the relation between the galaxy 
power spectrum Pgg{k) and the linear matter power spec- 
trum P\in{k)^ 



A^bin 


Mass Range 


C 


(3 




lO^^M© 


arcmin^/K 




Halo Bias 


781 


3.20 - 5.32 


1.41 ± 0.09 


1.51 ± 0.05 


470 


2.69 - 3.09 


1.48 ± 0.09 


1.31 ± 0.06 


1050 


1.76 - 2.05 


1.32 ± 0.05 


1.24 ± 0.04 


19003 


0.56- 1.04 


1.19 ± 0.03 


1.07 ± 0.02 


Weak Lensing Measurements 


3127 


1.72 - 3.46 


1.37 ± 0.03 


1.30 ± 0.02 


590 


1.62 - 2.48 


1.21 ± 0.07 


1.24 ± 0.05 


5408 


1.12 - 1.74 


1.16 ± 0.03 


1.20 ± 0.02 


11274 


0.60 - 1.02 


1.12 ± 0.03 


1.08 ± 0.02 


Note. 


— The two 


groups of mass ranges corres 


pond to the 



separate methods that are used to obtain halo masses for the LRG 
sample (see Table [T}. 



On the large scales probed by the LRGs, the bias b{k) 
is expected to be nearly scal e-independent, b ut str ongly 
dependent o n lum inosity. iTegmark et al.l (|2004f ) and 
IZehavi et al.l (|2005l ) have fitted an empirical model for 
galaxy bias relative to the bias of galaxiefl: 



b 



0.85 + 0.15— + 0.04(M^ - Mo.ir), (8) 



where - 51ogio(/i) = -20.44 (jBlanton et al.ll2003[ ). 
For h 0. 71, 

^•^r-band. iPercival et al.l (|2007 ^ have shown this model 
for relative bias to be a reasonable fit for the LRG sample 
considered in this pa per. The value f or used in this 
a nalysis is b^ = 1.19 (|R.eid et al.ll2010[ ). 

iTinker et al.l (|2QTQ^ use numerical simulations to de- 
rive a relationship between halo mass and bias. They 
identify halos in their cosmological simulation using the 
spherical over density algorithm to define the halo mass: 



Ma = -7rRlp{z)A, 



(9) 



where p{z) is the mean matter density of the universe at 
redshift z. This analy sis uses A = 2 00 and denotes the 
halo mass by M2oop- iTinker et al.l (|2010) estimate the 
relation between halo bias and M2oop- 



b{iy) = 1-A- 



0.183i/ 



1.5 



2.4 



(10) 



where u represents the "peak height" of the density field 
and is given hy u = Sc/a{M200p)^ is the critical density 
for collapse, and a is the linear matter variance on the 
scale of the halo, i.e., R = (3M/47rp)^/^. We use z = 
to determine the matter variance a in order to compare 
b{u) to b/b^. In all calculations, we use Sc = 1.686. The 



PM=b{kfPun{k). 



(7) 



is defined with respect to the Schechter (1976) luminosity 
function for SDSS galaxies (see Blanton et al. 2003) 
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parameters in equation [TOl are: 

A = 1.0 + 0.24^ exp[-(4/^)'^], 
5 = 0.44^-0.88, 

C = 0.019 + 0.107^ + 0.19 exp[-(4/^)'^], 

where y = log^^ol^OO). 

Using equations [8l and [T0\ we obtain estimates for halo 
mass for each luminosity bin. These values are given in 
Table [TJ We determine mass from luminosity for each 
LRG in a given bin and report M200P as the mean of 
these values. The error in M200P is estimated by the 
standard deviation of this mean, within each bin. 

Stacked gravitational lensing measurements are an al- 
ternative method of estimating the ha lo mass-galaxy lu- 
minosity relation. iReves et al.l (|2008f ) sta ck a sample of 
clusters from the SDSS maxBCG catalog (Koester et al. 
[2007l bl) and use weak gravitational lensing measure- 
ments to investigate the relation between cluster mass 
and v arious optical t racers. Since LRGs trace massive 
halos (|Ho et al.ll2009l ). the maxBCG catalog should rep- 
resent a population simil ar to the LRG sample used in 
this work. Noting that Re ves et aTl (!20b8) bin their clus- 
ter sample by ^-^^r-band BCG luminosity, we compute 
the ^-^^r-band luminosity for each LRG and de termine 
halo masses using equation 15(c) of iReves et al.l (j2308). 
Table [U lists these mass values, as well as those obtained 
from analysis of halo bias. The range in these two mass 
estimates represents some of the systematic uncertainties 
in relating LRG luminosity to halo mass. 

6. SZ-MASS SCALING RELATION 

In this section, we present an analysis of the correlation 
between cluster mass and integrated SZ signal. Specif- 
ically, we compare I200P to the cluster mass M200P for 
the luminosity /mass bins in Ta ble [J For self-similar evo- 
lution, analytical theory (e.g.. iKomatsu fc Seljak|[2001f ) 
predicts a simple relationship between Y and cluster 
mass. Assuming virial and hydrostatic equilibrium, the 
cluster gas te mperature can be relate d to M and E{z) 
through (e.g.. lBrvan fc NormanllT998f ) 

T (X M2/^£;(z)2/^ (11) 
where for a flat ACDM cosmology 

E{z) = [nm{i + zf + nA]^^^. (12) 

Under assumptions of an isothermal ICM, 

Y2oop(xdf{z) M200P-T, (13) 

where dA{z) is the angular diameter distance to the clus- 
ter. Combining the above equations leads to the expec- 
tation for self-similar evolution 

Y200P oc d^' M%1-^ E{zf'\ (14) 

Anticipating this scaling relation, we have plotted 
y2QQpd\E~'^^^ against the halo masses M200P in Figure [6l 
We show the two separate halo mass estimates for a given 
bin to illustrate the uncertainties that exist in convert- 
ing LRG luminosity to halo mass. The results for both 
the full (Table [2]) and radio-quiet (Table [3]) LRG samples 
are shown in this figure. The expected Y — M relation 



as determined from numerical simulations is shown as a 
solid line and is given by: 

where the best fit parameters are a = 1.76 and 7 
= 5.74. These parameters are obtain ed from the 
microwave sky simulations described in iSehgal et al.l 
(j2010f ). where halos of mass M200P > 2.82 x 10^^ M© 
have been considered. Furthermore, the redshift range 
of the clusters is restricted to 0.15 < z < 0.50, and a 
cylinder definition of Y is used. Tables [2] and [3] give 
the values of Y2oopd\E~'^^^ and associated errors for 
the full and radio-quiet samples. These estimates are 
determined using a summation similar to equation 
[6l where the redshift of each LRG is used to calcu- 
late dA{z) and E{z). They correspond to the masses 
derived from considering luminosity-dependent halo bias. 

7. DISCUSSION 

By stacking LRGs, we have obtained estimates of the 
SZ signal as a function of mass for groups and clusters. 
The stacking analysis enabled the detection of an SZ sig- 
nal for clusters with masses around lO^^M©, and the 
amplitude of the signal is consistent with our theoretical 
estimates. Using masses derived from halo bias, the slope 
of the inferred Y — M relation is consistent with expected 
results. However, when using stacked weak gravitational 
lensing measurements to estimate the halo masses, the 
slope appears to be slightly steeper than the self-similar 
prediction. 

For both mass calibrations, the overall normalization 
of the Y — M relation appears lower than expectedQ 
As the SZ effect traces the thermalized cluster pressure, 
there exists the possibility that there are significant non- 
thermal pressure components supporting the halo mass. 
However, the normalization of the Y — M relation de- 
pends upon the method of obtaining halo mass from LRG 
luminosity, which is subject to large uncertainties. Ap- 
plying an analysis of halo occup a tion d istribution to the 
SDSS LRG sample. IZheng et al.l (|2QQ9l ) find a -16% log- 
normal scatter in the central LRG luminosity for halos of 
a given mass. This scatter applies only to central LRGs 
and does not include the —5% of galaxies that are satel- 
lites. Thus, the scatter in the LRG luminosity for a fixed 
halo mass is likely even larger. Furthermore, there exists 
uncertainties in the mass estimation from weak lensing 
measurements (e.g. jMandelbaum et al.ir2008l : lRozo et al.l 
2009) that contribute to the discrepancy in our mass cal- 
ibrations. 

We have detected radio source contamination of the 
SZ signal at 148 GHz for -10% of the LRGs. This sug- 
gests the possibility that infrared sources may also have 
a contamination effect on our measurement of the SZ sig- 
nal for the LRG sample. Mittal et al. (2009) found that 

^ Note that immediately following the submission of this 
work, the Planck collaboration reported results on SZ scaling 
relations for both optical- and X-ray-selected galaxy clusters 
(Planck Collaboration et al. 2011a, b). Its analysis of the optical 
scaling relations finds a similar deficit in the overall normalization, 
while the scaling relations for X-ray clusters are in agreement with 
expected results. 
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the likelihood for a cluster to host a radio-loud BCG 
is a strong function of the central cooling time, in the 
sense that all strong cool-core clusters contain a radio- 
loud BCG. Thus, our radio-quiet LRG sample may select 
against clusters/groups with a strong cool core. However, 
in his an alysi s of a l ar^e sample of clusters and groups at 
z < 0.1, iSunI ()2QQ9l ) did not find any radio- loud central 
galaxies in groups with strong cool cores. If this is true 
at higher redshift, our radio-quiet LRG sample would be 
representative of all groups and clusters. 

Another source of systematic error is the conversion of 
the measured decrement STo to the integrated SZ signal 
Y. Ou r approach relies on simulations from Sehgal et al.. 
()2QlQf ) to estimate the effect of the matched filter on 
smoothing the cluster signal. Future work will explore 
this uncertainty by using simulations with different mod- 
els of cluster gas physics. 

This work is a first application of a novel stacking 
method for searching for the SZ effect in clusters by ap- 
plying a full optical catalog to a CMB survey. Both data 
sets are rich and have multiple components, with pos- 
sible correlations between them, that we will be sorting 
out as the data and our understanding improves. Future 
ACT analysis of the SZ effect from galaxy clusters will 
include data from multiple seasons and millimeter bands 
and will include the new SDSS DR8 analysis. 
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Bin 1 Bin 1 Bin 1 




(a) Full LRG Sample (b) Radio-quiet LRG Sample (c) Random Sample 



Fig. 2. — The stacked filtered temperature maps for each luminosity bin for the (a) full LRG sample, (b) radio-quiet LRG sample and 
(c) a randomly selected sample. Each map is 1' x 1' and is in CMB temperature units (|-tK). The values for bTo and the number of maps 
in each stack are given in Tables [2] and [3] for (a) and (b), respectively. The radio-quiet sample excludes the 10% of the full LRG sample 
that has been defined as radio-loud (see Section [STT}. The stacked random maps serve as a null test of the data, as the SZ signal should 
be consistent with zero for these maps. The number of maps in each stack for the random sample is equal to that of the full sample. The 
solid circles in each submap mark the characteristic size of a cluster for a given luminosity /mass bin, in terms of ^200- The values for ^200 
have been determined from simulation. The dashed circles show the FWHM of the ACT beam function. The crosses mark the center of 
each map. See online article for color version. 
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(b) Radio-quiet LRG Sample 



Fig. 3. — The distribution of the mean of 6To for each luminosity bin in the full and radio-quiet LRG samples, determined from bootstrap 
resampling. The mean and standard deviation a obtained from the resampled distribution are shown for each bin. The distribution of 
the mean is well-approximated by a Gaussian function, which is shown for each panel. The bootstrap errors a are reported in Tables [2] and 
[3] and are used throughout this work as the error estimate for STq. 



12 



Hand et al. 










-20 




-40 






o 


-60 




-80 




-100 




0.2 0.4 0.6 0.8 

simulated Y200P (10~^ arcmin^ ) 




0.1 0.2 0.3 0.4 0.5 

simulated Y200P (10~^ arcmin^ ) 



0.25 0.3 




simulated Yonn 7, (10 ^ arcmin^ ) 



0.04 

- 200p ' 



0.08 



simulated Yonn 7, (10 ^ arcmin^ ) 



Fig. 4. — The redshift depend ence of the relation between STo and Y2oop for each luminosity /mass bin, as determined from simulations 
described in'Sehgal et al.' mW). The redshift ranges are represented as follows: z < 0.26 (squares), 0.26 < z < 0.37 (circles), and z > 0.37 
(triangles). The four subplots correspond to the luminosity bins used in the LRG stacking analysis. The mass range for each bin corresponds 
to the mass estimates determined from analysis of halo bias (Table [1}. 
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Fig. 5. — The correlation between 5To and Y200p for each luminosity /mass bin, as determined from simulations described in lSehgal et al.l 

(|201C1I ). The simulation values for I20OP of the i^^ cluster are fit to the model Yest,i = C STo^i . The four subplots correspond to 
the luminosity bins used in the LRG stacking analysis. The mass range for each bin corresponds to the mass estimates determined from 
analysis of halo bias (Table 
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Fig. 6. — SZ - Mass Correlation. 1200/9 estimates are plotted against cluster masses for each luminosity bin for both the full (Table [2]) and 
radio-quiet (Table [3]) LRG samples . The two sets of points shown use separate mass estimates for the four LRG luminosity bins. The black 
circles correspond to mass estimates derived from analysis of halo bias. The triangles (red in online version) correspond to mass estimates 
derived from weak gravitational lensing measurements (Reves et al. 2008). The range in these masses represents some of the systematic 
uncertainties in converting LRG luminosity to halo mass. The solid line shows the expected model for the Y — M relation (equation I15|) . 
as determined from the microwave sky simulations of Sehgal et ah (2010. ). 



